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ABSTRACT 

A FORTRAN IV computer code for the micromechanics, macromechanics, 
and laminate analysis of multilayered fiber composite structural compo- 
nents is described. The code can be used either individually or as a 
subroutine within a complex structural analysis /synthesis program. The 
inputs to the code are constituent materials properties, composite geom- 
etry, and loading conditions. The outputs are various properties for ply 
and composite; composite structural response, including bending-stretching 
coupling; and composite stress analysis, including comparisons with fail- 
ure criteria for combined stress. The code was used successfully in the 
analysis and structural synthesis of flat panels, in the buckling analysis 
of flat panels, in multilayered composite material failure studies, and 
lamination residual stresses analysis. 
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INTRODUCTION 

Multilayered filamentary composites are structural materials which 
have evolved in recent years. These composites are constructed by lami- 
nating several plies (layers). The plies, in turn, consist of high 
strength filaments which are embedded in a low strength, low density 
(for nonmetallic matrices) matrix material while the matrix is in its 
molten state. Subsequently, the system is processed at specified pres- 
sures and temperatures. The result is a very light (for nonmetallic 
matrix composites), high strength material. Its strength depends on the 
orientation of the filaments with respect to the direction of the maxi- 
mum anticipated stress. Coincidence of the filaments with the direction 
of maximum stress utilizes this material most efficiently and gives it 
its most attractive features, stiffness and strength to weight ratio. 

Multilayered filamentary composites as structural materials have 
exhibited great potential in space vehicles, deep submergence vessels, 
radomes , and other structures where the stiffness and strength to weight 
ratio, nonmagnetic, structural damping, elevated temperature, and environ- 
mental resistant properties are important. The underlying structural 
strength principle for these materials is that the filaments are the pri- 
mary load carrying members while the matrix provides the structural in- 
tegrity by serving as a load transferring medium, providing rigidity 
(keeping the filaments fixed in their position) , and protecting the fila- 
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ments from exposure to unfavorable environment. 

Though multilayered filamentary composites have many desirable design 
features, as a structural material they are at the dawn of their great 
potential. One of the factors that keeps these composite materials at the 
dawn of their great potential is the designer's lack of proven and readily 
available mathematical formalisms which could be used to select constitu- 
ent materials when the structural application of the composite is speci- 
fied. These types of mathematical formalisms constitute a multilevel 
analysis encompassing the following levels: micromechanics, macromechanics, 

combined-stress strength functions, laminate analysis, and delamination 
criteria. This type of multilevel analysis is quite extensive and cumber- 
some. It needs to be transformed into an efficient computer code to be 
effective and practical. 

A FORTRAN IV computerized multilevel analysis for multilayered fiber 
composites which has grown to maturity over the past several years is the 
subject of this paper. This multilevel analysis includes all the levels 
mentioned previously. It has been extensively used for analyzing several 
fiber composite systems. The computer code is described herein from the 
engineer 's /analyst's usage viewpoint. Therefore, the description is 
limited to input/output, results produced and application versatility. 

The complete documentation of the code is given in [1]. Several of the 
mathematical functions programed are described briefly as part of the code 
output . 

OBTAINING AND USING THE CODE 

The mechanics required to use this code for the analysis of multi- 
layered fiber composites are described in some detail. Here, it is assumed 
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that the user is interested in using the code as a tool only and that he 
has available to him a FORTRAN IV manual. 

The physical representation of the code is illustrated in Fig. 1. 

The geometry of the constituents, the ply, and the composite are defined 
in this figure. The required input properties, correlation coefficients, 
and computed properties are summarized in Fig. 1 in symbolic form. 

The physical arrangement of the code is illustrated in Fig. 2. The 
numbers given in each block of cards are for subsequent discussion and do 
not appear on the code. Four steps are required to use the code in the 
user's computer facility: 

(1) Obtain the code. 

(2) Make it operational in the user's computer facility. 

(3) Supply the input data. 

(4) Interpret the code output results. 

Obtain the Code 

The code can be obtained in cards or tape from [2] . If this is not 
convenient or possible, then the cards can be punched from the compiled 
listing appended in [1], 

Make It Operational 

Making the program operational requires the availability of a FORTRAN 
compiler in the user's computer facility, certain control cards at the be- 
ginning of the code, and the card that preceeds each subroutine. Consult 
your computer group about these items. The control cards present in the 
code are only for the Lewis IBM 7044/7094 direct couple system. Once the 
deck of cards has been assembled as is shown in Fig. 2 (except Input Data) 
with the proper control cards, the user is ready to compile the code in 
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his facility. The compilation will indicate whether any additional modi- 
fications are needed. Most modifications will be minor and will usually 
deal with certain logical and subroutine call statements peculiar to each 
compiler. Consult your computer group for these modifications. 

Supply the Input Data 

The physical arrangement of the input data cards is illustrated in 
Fig. 3. The numbers in the group of cards are for identification purposes 
in this description and do not appear on the cards. A sample input data 
sheet is illustrated in Table 1 for the Thornel-50/ epoxy composite system. 

Listings of input data for several composite systems are appended 
in [1]. These systems are shown graphically in Fig. 4. The input data 
for these systems can be punched from the listings, and the cards that 
need alterations for the specific problem can be modified accordingly. 

Input data for additional composite systems may be easily prepared. 
This is done by selecting a related system from those available and modi- 
fying those entries that need modification. 

After the input data have been properly assembled (as is shown in 
Fig. 3), it is placed in its physical position (Fig. 2), and the code is 
ready to be run for results. 

Input Ply Properties 

There could be cases when the user would prefer to supply some of 
his own ply properties instead of using the code to compute them. The user 
has to provide his own formats for these cases. The physical location for 
these statements is described in the section MAIN PROGRAM. A flow chart of 
the MAIN PROGRAM is shown in Fig. 5. The functions of the subroutines are 


listed in Table 2. 
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OUTPUT AND MATHEMATICAL RELATIONS 
The program output consists of printing out (1) the input data, 

(2) the composite three-dimension strain-stress and stress-strain rela- 
tions about the structural axes (see Fig. 6) , (3) the composite properties 
generated in array PC, (4) the composite constitutive equations about the 
structural axes, (5) the reduced bending and axial stiffness, 

(6) displacement-force relations, (7) the current load or displacement 
condition, and (8) the ply properties, stresses, and margin of safety gen- 
erated in array PL. 

The printout of the input data is preceded by its code name. For 
example, the first and second lines of printout are 
THORNEL-50 /EPOXY 
NL ,NPL ,NPC ,NFPE ,NLC 
8 71 54 1420 1 

The output of code generated data is preceded by title headings. The 
output of the composite three-dimensional strain-stress temperature rela- 
tions and composite stress-strain relations about the structural axes are 
printed under the headings 

3-D COMPOSITE STRAIN STRESS RELATIONS - STRUCTURAL AXES 

The matrices [E ] ^ and {a } in the equation 
c s c s M 

(e } = [E ] _1 {a } - AT {a } 

c s c s c c s 

are printed out. 

3-D STRESS STRAIN RELATIONS - STRUCTURAL AXES 

The matrix [E ] in 
c s 

(a } = [E ] {e } 

c s c s c s 

is printed out. The subscript s in the preceding equations indicates 
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that the relations are written about the structural axes. It is noted 
that these properties are only local to subroutine GACD3 . They can be 
made global if needed. 

The output of the composite properties, generated in array PC are 
printed under the heading 

COMPOSITE PROPERTIES - VALID ONLY FOR CONSTANT TEMPERATURE THROUGH 
THICKNESS 

LINES 1 TO 31 3-D COMPOSITE PROPERTIES ABOUT MATERIAL AXES 
LINES 33 TO 54 2-D COMPOSITE PROPERTIES ABOUT STRUCTURAL AXES 
Fifty-four entries are printed under this heading as follows: 


PC(1) 


P c 

weight density 

PC (2) 


t 

c 

thickness 

PC(3) 

to PC(ll) 

ny 

three-dimensional stress-strain re- 
lations about material axes 

PC(12) 

to PC (14) 

(a } 
c 

three-dimensional coefficients of 
expansion about material axes 

PC (15) 

to PC (18) 

{ V> H c 

three-dimensional heat conductivities 
and heat capacity along material axes 

PC(19) 

to PC(30) 

E cll’ G cl2 ,V cl2 

three-dimensional constants about 
material axes 

PC(31) 


z 

distance to reference plane from bottom 
of composite 

PC(32) 



blank 

PC (33) 

to PC (38) 

t E c ] 

two-dimensional stress-strain relations 


about structural axes 
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PC(39) to PC(47) E ,.,G -,~,v , „ two-dimensional elastic constants 

ell’ cl2 cl2 

along structural axes 

PC(48) to PC(54) {a c ),K c ,H c two-dimensional coefficients of 

thermal expansion, heat conductiv- 
ities, and heat capacity along 
structural axes 

The output for the composite constitutive equations are printed under 
the heading 

FORCES FORCE DISPLACEMENT RELATIONS DISPL THERMAL FORCES 



The elements of matrices A , C , N , and M are printed out. 

cx ex’ cATx cATx K 

The output for the reduced bending rigidities is printed under the 
heading 

REDUCED BENDING RIGIDITIES 

R 

The elements of D are printed out in one line. 

cx 

The output for the reduced axial stiffness A is printed out under 

CX 

the heading 

REDUCED STIFFNESS MATRIX 

The inverse of the constitutive equations is printed out under the 
heading 


DISP DISPLACEMENT FORCE RELATIONS FORCES 
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The elements of this inverse are printed out. 

The output for the current load condition is printed next to the 
headings 

FOR THIS CASE NBS(X,Y,XY-M) IS 
and 

FOR THIS CASE MBS (X ,Y ,XY-M) IS 

The current values of N , N , N , M , M , and M are printed 

cx cy cxy cx cy cxy 

out under these headings. 

The output for the current displacement conditions is printed under 
the heading 

FOR THIS CASE THE DISPLACEMENTS DISV (ECSXX ,ECSYY ,ECSXY ,WCBXX,WCBYY ,WCBXY) ARE 
The output of the ply properties generated in array PL are printed out 
under the heading 

LAYER PROPERTIES, ROWS -PROPERTY , COLUMNS-LAYER 

Seventy-one entries are printed out under this heading as follows: 

PL (1,1) ply void content 

PL(2,I) kf ply apparent fiber content 

PL(3,I) k^ ply actual fiber content 

PL (4, I) k^ ply apparent matrix content 

PL (5, I) k^ ply actual matrix content 

PL (6, I) ply weight density 

PL(7,I) t^ ply layer thickness 

PL (8, 1) <5 ply and interply layer thickness 

PL(9,I) H_. interply layer distortion energy coefficient 

PL(10,I) "z distance from bottom of composite to ply centroid 

z distance from reference plane to ply centroid 


PL(11,I) 
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PL(12,I) 0 angle from structural axes to 

cs 

composite material axes (same 
for all plies) , Fig. 6 

PL (13, I) 0 angle from ply material axes to 

X/ 

composite material axes, Fig. 6 

PL (14, I) 8 angle from ply material axes to 

x»S 

composite structural axes, Fig. 6 


PL(15 ,1) 

to 

PL (23 ,1) 

'V 

ply stress-strain relations 

PL(24,I) 

to 

PL(26,I) 

<»t> 

ply thermal coefficients of expansion 

PL (27 ,1) 

to 

PL(29 ,1) 

<V 

ply heat conductivities 

PL (30, I) 



H cjl 

ply heat capacity 

PL(31,I) 

to 

PL(32,I) 

E m ,v n2 ,G £i2 

ply elastic constants 

PL(43,I) 

to 

PL (48, 1) 

P y22 ,P yl2 ’ p yl3 

ply strain magnification factors 

PL(49,I) 



P yde£ 

interply delamination factor 

PL(50,I) 



AT 

ply temperature 

PL(51,I) 

to 

PL(60 ,1) 

S £11T’ StC ' 

ply limiting stresses 

PL (61 , I) 



K £12ag 

coefficient in combined-stress - 





strength criterion 

PL (6 2, I) 




ply margin of safety — from dis- 





tortion energy 

PL (63, I) 




interply delamination margin of 





safety 

PL(64,I) 

to 

PL(69 ,1) 

{£ * } - ta t } 

ply applied strains and stresses 

PL (70 , I) 



A P . 

adjacent ply relative rotation 

PL(71 ,1) 




margin of safety from Hoffman's 


failure criterion 

The failure criterion may be determined by either of the following methods. 
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(1) Modified distortion energy 


F = 1 - 


(r 1 ^) 

V 3 2,11a/ 


'me \ 


£223 


- K 


,K, 


°£lla °£22 


£12a3 112«B S mo S, 22 


£12S A 


£12S 


1 


PL(62,I) 


-*1 


(2) Hoffman's criterion 

2 


F = 1 - 


°m °£ii a £22 


422 


S mC S 411T 


S £22C S £22T 


, S 411C S U1T 

s s °jni 

me 411T 


S £22C S £22T °£12 

S „_ S _ °£22 + „2 


£22C £22T 


£12SJ. 


PL(71,I) 


F > 0 no failure 
F = 0 incipient failure 

F < 0 failure 

where a denotes stress and S the limit stress always taken positive. 
The subscripts denote direction and sense. 

th 

The interply delamination criterion for the j interply layer at the 
til 

m load condition is governed by 

-*■ PL (63, 1) when i > 1 


! - (ML 

Wdel 


JJ 

A<j>. = -r (e - e ) (sin 20. - sin 20. ,) + -q- e (cos 20. - cos 20. ,) 
Y j 2 cyy exx l i-1' 2 cxy v l i-1 7 


(e } = [A J _1 ({ N } + {N A _ } + [C ]{w , }) 
cx cx \ cx cATx 1 cx cbx / 


or as given by the displacement force equation described previously. 

The inputs to the subroutines to generate the PL array are the 
ply angle measured from the structural axes (0 from PL(14,I)), 
the distance from the reference plane 
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to centroid of the ply (z^ from PL(11,I)), the ply temperature (AT^ 

from PL(50,I)), the interply delamination limit (A$^ e ^ from PL(60,I)), 

and the ply thermoelastic properties stored in PL (24 to 26,-1) and 

PL(31 to 42,1). The ply extensional and coupling rigidities A = ACX 

cx 

and C = CPC; the local curvatures w , = WXX; the adjustment constants 

cx cbx J 

K Il2TT * BETC1 - 7 >- K U2CT ‘ BEI(2 ' 7 >- K U2TC ‘ “ d 

k' = BET(2,8); and the load conditions N = NBS(m). All other ad- 

cx 

justment constants are input data. 

A sample output for a Thornel-50/ epoxy angle ply composite is given 
in Table 3 . 

CODE APPLICATIONS 

This computerized multilevel analysis has been applied extensively 
to various aspects of fiber composite analysis, design, and optimization. 
Typical cases include: nonmetallic [3,4,5] and metallic [6] fiber com- 

posite characterization, random composite characterization [7] , lamination 
residual stresses [8,9,10], buckling of fiber composite anisotropic plates 
[11,12], exploratory designs of fiber composite turbine blades (unpublished 
notes) , feasibility studies of high-tip-speed compressor blades made from 
advanced fiber composites [10] , failure analysis of flat and tubular 
members [3,4], gross and local impact studies [13,14], structural synthe- 
sis, and sensitivity studies of flat plates [3], and other important fac- 
tors in designing with fiber composites [15] . 

IMMEDIATE EXTENSIONS 

The code can be readily modified and supplemented to handle nonlinear 
material response, temperature dependent properties, inter, and intraply hybrid 
fiber composites. The details of these modifications become apparent once 
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the user has some experience in using the code. 
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TABLE I. - INPUT DATA FOR THORNEL- 50/E POXY COMPOSITE 


THORNEL-50/ EPOXY 


8 71 54 

1420 1 




.50000E.08 

.10000E.07 

. 1 OOOOE . 07 

. 20000E .00 

•25000E.OO 

.20000E .00 

. 1 3000E .07 

. 70000E .06 

. 130UOE.07 

• 57000E • 06 

• 57000E • 06 

•57000E.06 

• 36000E • 00 

. 36000E .00 

. 360OOE.O0 

•00000E.00 

.OOOOOE.OO 

•OOOOOE.OO 



.40000E.01 

•20000E.01 

. 40000E . 0 1 

.20000E.01 

.OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•23560E.01 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

. 10000E .01 

•10000E.01 

. 10000E.01 

•OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

-.55000E-06 

•56000E-05 

•56000E-05 



•42800E-04 

•42800E-04 

•42800E— 04 



•58000E.03 

•58000E.02 

• 58000E . 02 

. 17000t.00 

.12500E.01 

• 12500E .01 

. 12500E.01 

• 25000E . 00 

•OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

•22300E .00 




•10000E.01 

. 10000E.01 

. 10500E.01 

.10500E.01 


•31416E.01 

F 

F 





F 

F 

.OOOOOE.OO 

. 05900E • 00 

• 04430E .00 

.00026E.00 


. OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

O.OOOOOE.OO 

O.OOOOOE.OO 

0» OOOOOE .00 



0.50000E-00 

O.50OOCE.O0 

0. 50000E • 00 

0.50000E.0O 

0.50000E.OO 

0 . 50000E .00 

0.50000E.00 

0.50000E .00 



•OOOOOE.OO 

•45000E.02 

-. 45000E .02 

•90000E.02 

•90000E.02 

-.45000E.02 

•45000E.02 

•OOOOOE.OO 



•0O8O5E.OO 

.00805E.00 

. 00800E . 00 

•U08O0E.00 

.OOBOOE.OO 

•00805E.OO 

.OO8O5E.O0 

• 00800E .00 



-0.30000E.03 

-0.30000E.03 

-0. 30000E .03 

-Q.30000E.03 

-0.30000E.03 

-0.30000E.03 

-0.30000E.03 

-0. 30000E .03 



•83000E.00 

. 10000E.01 

• 26000E .00 

. 2 7000E . 00 

.17OOOE.O0 

• 16500E .02 

. 10000E.01 

. 10000E.01 

• 04650E .00 

. 10000E .01 

.50000E.O0 

.13300E.02 

•31900E.05 

. 10000E .01 

. 10000E.01 

. 1 OOOOE • 0 1 





.23000E.06 

•21000E.05 

• 02000E . 00 

.050OOE.0O 

. 04500E . 00 

•04500E.00 





. 50000E .04 

•OOOOOE.OO 

•OOOOOE.OO 



. 10000E.03 

.OOOOOE.OO 

•OOOOOE.OO 



O.OOOOOE.OO 

O.OOOOOE.OO 

O.OOOOOE.OO 

O.OOOOOE.OO 

O.OOOOOE.OO 


O.OOOOOE.OO 



TABLE 2. - SUBROUTINE DESCRIPTION 


Subroutine 


Function 


INVA 

GLLSC 

GACD3 

GPCFD2 

GPH 

GECL 

GSMF 

COMPS A 


Matrix inverse 

Simple limit stresses (strengths) ply (unidirec- 
tional composite) 

Three-dimensional laminate thermoelastic properties 
Two-dimensional laminate thermoelastic properties 
Ply heat conductivities 
Ply thermoelastic properties 
Strain magnification factors 

Laminate analysis, ply stress-strain, ply margin of 


safety 



TABLE 3. - SAMPLE CASE OUTPUT 

Angle Ply Composite 


TH0RNa-5C/CPOAY 


NL.NPUNPC.NFP6 »NLC 

6 7 1 5 * 1**0 1 


fcf 11. EF22#fK 23, NUM2 .NUF 23 ,NUF i 3 ,EF 1* ,6F23 ,EF13»6M11 , 6H2 2 , EM33, NUM 12. NU* 23 NUM13.EM 12, EH ’5,6 HI 3 


0.500001 OH 

'j. no oc 6 

07 

3. 1 00006 

07 

0.200036 

33 

0. 250006 

03 

0.20000E 

00 

0.130006 

07 

0. 70000= 

06 

0.13300E 07 

3.573036 36 

0.57CQ0E Uo 
VCf 

0. 570006 

06 

0.360006 

00 

0.360006 

30 

0.36000E 

03 

0. 


0. 


c. 




o.*cccot 01 

0 • 200306 

01 

0.400006 

01 

0.200006 

01 

0. 


0. 


0. 


0.235606 

01 

0. 

3. 

O'. 

0. 


0. 


C. 


O.IOOOOE 

01 

0. 10000E 

01 

0. 10000E 

01 

0. 


0. 

3. 

VAF 

-0.55CCOE-C6 
V AM 

f). 560 006- 

OS 

0. 5600-36 

-05 













0. *28006-0* 
CH* 

0.428006- 

0* 

0.428006 

-0* 













0°. 560006 0 1 

0.5800C6 

02 

0.5 30C06 

02 

0.170006 

00 

0.125006 

01 

UJ 

o 

o 

IT 

CM 

o 

01 

0.125006 

01 

0.250006 

00 

0. 

0. 


0. 0.225CC6 00 

6T A 

0.10000E 01 0.100006 01 0.105006 01 0.10S00E 31 

pie 

0.3l*16fc 01 

TL INP 
F 

CSANd 

F 

b ICE 

F 

ft INOV 
F 


T HCS .RHflP.R HCH.UUF 


0. 


•3.53O0CE- 

01 

0.4*3006 

-01 

0. 260006-33 



KVL 










0. 


0. 


0. 


0. 


0. 


KFL 










0.5000 06 

0 7 

0.5000CE 

00 

0.50000E 

00 

0. 50000E 

03 

0.50000E 

03 

TMLC 

0.3CC006 

TL 

02 

-J. 2-30CCE 

02 

0.30000E 

02 

UJ 

o 

o 

o 

o 

0 

1 

02 

-0.30 000 E 

02 

0 . 80 500E-02 

3 . 8050C6- 

02 

0. 80000E 

-02 

0. 80030 E -3 ?. 

0.800006-02 

PTEMP 

-0.300006 

0 3 

- ). 30000E 

03 

-0. 300 OOE 

03 

-0; 300006 

33 

-0.300006 

03 

BET 










C.82COOE 

CC 

0. lUJOCC 

01 

C. 2600CE 

00 

0. 27003 E 

33 

0.1 7000 E 

03 

0.50000E 

00 

0. 13300E 

02 

0. 31 9006 

05 

0. 13333E 

31 

0. 13033 E 

01 

L SC 










C.2300CE 

oe 

U. 210006 

05 

0. 200006 

-01 

0. 50000E-01 

0. *50006-01 


0. 0. 0. 

0.500006 00 0.50000E CO C. 500006 00 

0.30000E 02 - 0 • 300C0E 02 0.30000E 02 

0. 805006-02 C . 80500E- 02 0. 800006-02 

-0. 30000E 03 -0. 300006 ^3 -0. 30000E 03 

0.165306 02 0 ■ 100 COE 01 0. 100006 01 0.46500E-DI 3.133036 01 

0. 1003 OE 01 

0. *50006-01 


N8S 

0.50000E 0* 0. 0. 


MBS 

0,500006 02 0. 0. 

0I$V1 

0, 0. 0. O. 0. 0. 



TABLE 3. - SAMPLE CASE OUTPUT (CONT'D) 

3-D COMPOSITE STRAIN STRESS TEMPERATURE RELATIONS - STRUCTURAL AXES 


0. LA3UE-06 

-0.2449E-06 

0.4*716-37 

0. 

0. 

0. 

0.2449E-U6 

0. 7909E -Ob 

-0. 22466-06 

- 0 . 

-0. 

-0. 

0.4471E-07 

-0.2246E-06 

0*93536-06 

0* 

0. 

0. 

0. 

0. 

0 . 

0 .2937E-05 

0 . 

c. 

0. 

0. 

0 . 

0 . 

0. 2937E-05 

0. 

0. 

0 . 

0 . 

0 . 

0. 

0.2000E-06 


3-D COMPOSITE STRESS STRAIN RELATIONS - STRUCTURAL AXES 


-0. 31 31E-05 
0. II99E-0* 
0.2933E-04 
- 0 , 

- 0 . 

- 0 . 


0.15 04E 

OB 

0. 4 7 79 E 

07 

0.4285E 

06 

-0. 

-0. 

0. 

0.47796 

07 

0.28756 

07 

0.4620E 

06 

-0. 

-0. 

0. 

0.42856 

06 

0.4620E 

06 

0. 1160E 

07 

-0. 

-0. 

-0. 

0. 


-0. 


-0. 


0.3405E 06 

-0. 

-0. 

0. 


-0. 


-0. 


-0. 

0.3405E 06 

-0. 

3. 


0. 


-0. 


-0. 

-0. 

0.5001E 


COMPOSITE PROPERTIES - VAL/O ONLY FOR CONSTANT TEMPERATURE THROUGH THICKNESS 
LINES I TO H 3-0 CUMPUSITE PROPERTIES ABOUT MATERIAL AXES 
LINES 33 TU 5* 2-U COMPUSI TE PROPERTIES ABOUT STRUCTURAL AXES 


1 

RHOC 

0.51 65E- 

'01 

2 

TC 

0.640CE- 

-01 

3 

ecu 

0.1504E 

08 

4 

CC 12 

C.4779F 

07 

5 

CC 1 3 

C.4285E 

06 

6 

CC 22 

0.28756 

07 

7 

CC 23 

C.462CE 

06 

8 

CC 3 3 

C.U60E 

07 

9 

CC44 

C. 34056 

06 

10 

CC 55 

0.3405E 

06 

n 

CC 66 

C • 500 IE 

07 

12 

CT61 1 

-0.31 31E-05 

13 

CTE22 

0.11 99F-04 

14 

CT63J 

0.29 33E- 

'04 

15 

HU 1 

0.21B96 

03 

16 

HC 22 

0. 7644F 

02 

17 

HR 33 

0.3715E 

01 

18 

PMC 

0.2U43E 

00 

19 

ECU 

0.6992E 

07 

20 

EC 22 

C.U64E 

07 

21 

EC 33 

0.1069E 

07 

22 

EC 23 

C.3405E 

06 

23 

EC 31 

0.3405F 

06 

26 

ECU 

C. 50016 

07 

25 

NUC12 

0.1712E 

01 

26 

NUC21 

0.3096E 

00 

27 

NUC13 

- 0.3126E 

00 

28 

NUC31 

-0.4790E-01 

29 

NUC2 3 

0.2839E 

CO 

30 

NUC32 

0, 24 OIF. 

00 

31 

2CGC 

0.3200E- 

■01 

32 

B20EC 

0. 


33 

ecu 

o.i49eE 

08 

34 

ecu 

0.46006 

07 

35 

CC 13 

0. 


36 

CC 22 

0.2691E 

07 

37 

CC23 

C. 


38 

CC 33 

C . 500 IE 

07 

39 

ECU 

G.6992E 

07 

40 

EC 22 

0.12646 

07 

41 

EC 12 

0.5001E 

07 

42 

NUC12 

0.1712E 

01 

43 

NUC21 

0.30966 

00 

44 

CSN13 

-C. 


45 

CSN3I 

-0. 


46 

CSN23 

-0. 


47 

CSN32 

0. 


48 

CT£1 l 

-0.31 31E- 

05 

49 

CTE22 

0.U99E- 

04 

50 

CTE12 

-0. 


51 

MK 1 1 

C.2189E 

03 

52 

ML 22 

0.7544E 

02 

53 

ML 12 

-0. 


54 

HHC 

0.2043E 

00 



TABLE 3 


SAMPLE CASE OUTPUT (CONT’D) 


FORCES 



FORCE DISPLACEMENT RELATIONS 



01 SPL 

THERMAL FORCES 

NX 

0.9525E 06 

0.29*96 06 

0. 

0. 12216-03 

0.45786-0* 

0. 1068E-03 

ux 

-0. 1666 E 03 

NY 

0.29*96 06 

0.1722E 06 

0. 

0 .*578E-0* 

0. 2289E-0* 

0. 38156-04 

vx 

-0.3*276 03 

NXY 

0. 

0. 

0.32006 06 

0.10686-03 

0. 38156-04 

0. 30526-04 

VXPUY 

0. 

MX 

0.12216-03 

0. *5786-0* 

0.1068E-03 

0.3251E 03 

0.10076 03 

0.63906 02 

WXX 

-0. 22356-07 

MY 

0. *5786-0* 

0.2289E-0* 

0. 3Q15E-0* 

0. 10076 03 

0.58796 02 

0.2260E 02 

WYV 

-0.2980E-07 

MXY 

0 . 10686-03 

0. 3815E-0* 

0.30526-0* 

0 .6 393 E 02 

0.22606 02 

0.1092E 03 

MXY 

0.29806-07 


REDUCED BENDING REG IOI TIES 

0.225136 03 0.10067E 03 0.639006 02 Q.10067E 03 0.58793E 02 0.22595602 0.63900E 02 0.225956 02 0.1092*6 03 


RECUCEO STIFFNESS MATRIX 

0.9*2526 06 0.29*9*6 06 -0. *92336-10 0.29*9*6 06 0.172256 06 -0.192816-10 -0.492336-10 -0.192816-10 0.320036 06 


01SP. 


DISPLACEMENT FORCE RELATIONS 


FORCES 


UX 

0.22 3 56 - 05 

-0. 3827E-05 

0.11326-21 

-0.3*566-12 

0.64096-12 

-0.7792E-12 

NX 

vx 

-0.38276-05 

0.12366-0* 

0.15596-21 

0. 62256-12 

-0.27566-11 

-0.36806-12 

NY 

VXP8Y 

0.11326-21 

0.15596-21 

0.31256-05 

-0 .8006 E-12 

-0.54436-12 

-0. 2920E- 12 

NXY 

MXX 

- 0.3*566-1 2 

0.62256-12 

-0.80066-12 

0. 68306-02 

-0.11046-01 

-0. 17126-02 

MX 

WYY 

0.6*096-12 

-0.2756E-11 

-0.5**36-12 

-0.11046-01 

0.36316-01 

-0.1055E-02 

MY 

MXY 

-0.7792E-1 2 

-0. 3680E-12 

-0.29206-12 

-0.1712E-02 

-0 . 105 5E-C2 

0.10376-01 

MXY 


DISP. DISPLACEMENT FORCE RELAT IONS FORCES 


0.12116-01 

0.22356-05 

-0.38276-05 

0.U326-21 

-0.3*566-12 

0.64096-12 

-0.77926-12 

0.48336 

0* 

-0.22736-01 

-0.382 76-05 

0.12366-0* 

0.15596-21 

0.62256-12 

-0.27566- l 1 

-0.36806-12 

-0.3*276 

03 

-0.40036- 10 

0.11326-21 

0.15596-21 

0.31256-05 

-0. 8006E- 12 

-0.5*436-12 

-0.29206-12 

0. 


0.3415E 00 

-0.3*566-12 

0.62256-12 

-0.80066-12 

0.6833E-02 

-0.1 1046-01 

-0.17126-02 

0.50006 

02 

-0.5518E 00 

0.6*096-12 

-0.27566-11 

-0.5**36-12 

-0.110*6-01 

0.36316-01 

-0. 10556-02 

-0.29806 

-07 

-0.85616-01 

-0.77926-12 

-0.36806-12 

-0.29206-1 2 

-0.17126-02 

-0.1C55E-C2 

0.103 76 - 01 

0.29806 

-07 


FOR 

THIS 

CASE 

N 8 SI X * V » XY-M) 

1 S 

5000. 

0 . 

0 . 

FOR 

THIS 

CASE 

MBS(X'V.XY-M) 

1 S 

50. 

0 . 

0 . 


LAYER PROPER HES. ROKS-PRGPERTY » COLUMNS-L AYER 


1 

KV 

0. 

0. 

0. 

0 . 

0. 

0. 

0. 

0. 

2 

KF 

0. 50 006 00 

0.50006 00 

0.50706 00 

0.50306 30 

0.50006 00 

0.500C6 00 

0.50006 33 

0.53306 77 

3 

KFB 

0. 50 OOF 00 

0.5000E 00 

0.50036 00 

0.50306 00 

0 . 50 COE 00 

0.50006 00 

0.50006 37 

0.53006 03 

* 

KM 

0.50006 00 

0.50006 00 

3.50006 00 

0 . 533 36 30 

0.50 COE 00 

0.5000E "0 

0.50006 03 

0.53036 73 

5 

KM 8 

0.50006 00 

0.50006 00 

0.50336 00 

0.500C6 00 

0.50006 00 

0. 50036 00 

0.50006 00 

0.5300E 33 

6 

RMOL 

0. 51656-01 

0.5165E-01 

3. 51656-01 

0 .51656-01 

0.51656-01 

0.51656-01 

0.51656-01 

0.51656-31 

7 

tl 

0. 8O00E - 02 

0. 80006-02 

0.83336-02 

0.80006-02 

0.80C06-02 

0.90006-02 

0. 800 06-02 

0. 97006-32 

8 

OELTA 

0.658o6-0* 

0.65866-0* 

0.6 586 6-0* 

0.65866-0* 

0.65866-0* 

0.65066-0* 

0.65966-0* 

0.55956-3* 

9 

1L0C 

-0.00006-19 

0. 

0. 

0 . 

0 . 

0. 

n. 

0 . 

10 

ZB 

0.40006-02 

0. 12006-01 

0.2003E-01 

0 • 283 0 E-0 1 

0.36006-01 

0.44006-01 

0. 52006-01 

0.53006-31 

11 

ZGC 

-0.28006-01 

-0.20006-01 

-3.12036-31 

-0 .40306-02 

0.40006-02 

0.1 2006-01 

0. 23006-31 

0.29036-01 

12 

TMCS 

0. 

0. 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

13 

TMLC 

0.52366 00 

-0.5236E 00 

0.52366 OO 

-0 . 523 66 00 

-0.5236E 00 

0.52366 00 

-0.52366 30 

0.5236E 37 

l* 

THLS 

0.52366 00 

-0.52366 00 

0.52366 00 

-0.5236E 00 

-0.52366 00 

0.52366 00 

-0.52366 30 

0.52356 33 

15 

SC11 

0.25*96 08 

0.25*96 08 

0.25*36 0? 

0.25*96 08 

0.25*96 08 

0.25*96 08 

0.25*96 09 

0.25*96 38 

16 

SC 12 

0.41186 06 

0.41186 06 

0.41186 06 

0.41186 36 

0.4U8E 06 

0.4118E 06 

0.41186 06 

0.41186 06 

17 

SC 13 

0.41186 06 

0.41186 06 

0.*1U6 06 

0.41186 06 

0.4U86 06 

0.41186 06 

0.41186 05 

0.41186 05 

18 

SC 22 

0. 11606 07 

0.11606 07 

0.11676 07 

0.11606 07 

0.1160E 07 

0.11606 07 

0.11636 07 

O.1160F 37 

19 

SC 23 

0.47876 06 

0.47876 06 

0.4737= Ofc 

3 .478 7E 06 

0.4787E 06 

0.4 7 0 76 0 6 

0.47876 05 

0.47876 06 



TABLE 3 


SAMPLE CASE OUTPUT (CONC'D) 


20 

SC 33 

0.H60E 07 

0.1160E 07 

0.11636 07 

0.11606 07 

0.1160E 07 

0.1160E 

07 

0.11606 

07 

0.11606 07 

21 

SC44 

0.3405E 06 

0. 3405E 06 

0.34056 06 

0.34056 06 

0.3405E 06 

0.3405E 

06 

0.3405E 

06 

0.34056 06 

22 

SC 55 

0.34056 06 

0.3405E 06 

0. 3405 E 06 

0.34056 06 

0.3405E 06 

0.3405E 

06 

0.34056 

06 

0.34056 06 

23 

SC66 

0.6339E 06 

0.63396 06 

0.63396 06 

0.6339E 06 

0.63396 06 

0.63396 

06 

0.63396 

06 

0.63396 06 

24 

CTEll 

-0.6138E-07 

-0.61 38E-07 

-0. 61386-07 

-0.61386-07 

-0.61386-07 

-0.6138E- 

-07 

-0.61386- 

-07 

-0.61386-07 

25 

CTE22 

0.233 4E - 04 

0.23346-04 

0.23346-04 

0.23346-04 

0.23346-04 

0.23346- 

-04 

0.23346- 

-04 

0.23346-04 

26 

CTE33 

0. 2334E-04 

0.2334E-04 

0. 2334E-04 

0.23346-04 

0.2334E-04 

0.2334E* 

-04 

0.23346- 

-04 

0.23346-04 

27 

HKU 

0. 2906E 03 

0.29066 03 

0.29066 03 

0 . 290 66 03 

0.2906E 03 

0.2906E 

03 

0.29066 

03 

0.29066 03 

28 

HC 22 

0.37156 01 

0.37156 01 

0.37156 01 

0.37156 01 

0.3715E 01 

0.371 5E 

01 

0.3715E 

01 

0.3715E 01 

29 

H* 33 

0.37156 01 

0.37156 01 

0.3715E 01 

0.3715E 01 

0.3715E 01 

0.3715E 

01 

0.3715E 

01 

0.3715E 01 

JO 

KL 

0.2043E 00 

0.20436 00 

0. 2043 E 00 

0.2043E 00 

0.2043E 00 

0.2043E 

00 

0.20436 

00 

0.20436 00 

31 

EL 11 

0.25286 08 

0.2528E 08 

0.2528 E 08 

0.2528E 00 

0 .2 5? 8E C8 

0.25286 

08 

0.252 86 

08 

0.25286 08 

32 

EL 22 

0.959TE 06 

0.95976 06 

0.9597 E 06 

0.95976 06 

0.95976 06 

0.95976 

06 

0.95976 

06 

0.35976 36 

33 

El 33 

0.9597E 06 

0.95976 06 

0.95976 06 

0.9597E 06 

0.95976 06 

0.9597E 

06 

0.9597E 

06 

0.95976 06 

34 

GL 23 

0.J4C56 06 

0.34056 06 

0.34056 06 

0 .3435E 06 

0.3405E 06 

0.340 56 

06 

0.34056 

06 

0.34056 06 

35 

Gl 13 

0.6339E 06 

0.63396 06 

0. 6 339 E 06 

0 .63396 06 

0 .63396 06 

0.63396 

06 

0.63396 

06 

0.6339E 06 

36 

a 12 

0.6339E 06 

0.63396 06 

0.63396 06 

0 .63396 06 

0 .6 3396 06 

0.63396 

06 

0.6339E 

06 

0.53396 06 

37 

NUl 12 

0.2514E 00 

0.2514E 00 

0.2514E 00 

0.251 4E 00 

0.25146 00 

0.251 46 

00 

0.2514E 

00 

0.25146 00 

38 

NUL21 

0.9541E-02 

0. 95416-02 

0.95416-02 

0.9541E-02 

0 .954 IE- 02 

0.95416 

-02 

0. 95416 

-02 

0.95416-02 

39 

NUL13 

0.2514E 00 

0.25146 00 

0.25146 00 

0.2514E 00 

0.25146 00 

0.251 46 

00 

0. 25 1 46 

00 

0.25146 00 

40 

NUl 31 

0.95416-02 

0. 9541E-02 

0.95416-02 

0.954 IE-02 

0.95416-02 

0. 954 IE 

-02 

0. 9541E 

-02 

0.9541E-02 

41 

NUl 2 3 

0.4094E 00 

0.4094E 00 

0.40946 03 

0.4094E 00 

0.40946 00 

0.40946 

00 

0.40946 

00 

0.40946 00 

42 

NUL 3 2 

0.4094E 00 

0.4094E 00 

0.40946 00 

0.4094E 00 

0.4094E 00 

0.40946 

00 

0.40946 

00 

0.40946 00 

43 

SHFK22 

0.1522E 01 

0.15226 01 

0.15226 01 

0.152 2E 01 

0.1522E 01 

0.15226 

01 

0.1522E 

01 

0.15226 01 

44 

SMFD22 

0. 1918E 01 

0.19186 01 

0.19186 01 

0.1918E 01 

0.19186 01 

0. 19186 

01 

0.1 9186 

01 

0.1918F 01 

45 

SHFS 22 

0.1344E 01 

0.13426 01 

0.14086 01 

0.14306 01 

0.13996 01 

0.13976 

01 

0. 13996 

01 

0.1395E 01 

46 

SHFC22 

- 0. 0003E - 1 9 

-0.00006-19 

-0.33336-19 

-0.0030E-19 

-0.00006-19 

-0.00006 

-19 

-O.OOOfE 

-19 

-0.00006-19 

47 

SHF S 12 

0.3C24E 01 

0.30246 01 

0.33246 01 

0.30246 01 

0.30246 01 

0.302 46 

01 

0. 302 46 

01 

0.3024E 01 

48 

SHFS23 

U. 1396E 01 

0 • 1 3 96 6 01 

0.13966 01 

0.13966 01 

0.1396E 01 

0. 13966 

01 

0.13966 

01 

0. 1 3966 01 

49 

1LMFC 

-0. OOCOE-l 9 

0. 70606 02 

0. 7063 E 02 

0 .7060E 02 

0.70606 02 

0.70606 

02 

0.70606 

02 

0.706.06 02 

50 

TEMPO 

-0.3000G 03 

-0.30006 03 

-0. 3333 E 03 

-0.30006 03 

-0.30006 03 

-0.30006 

03 

-0.300CE 

03 

-0.30006 03 

51 

L SCI IT 

0.96 76E 05 

0. 96 766 05 

0.96766 05 

0.96766 05 

0.96766 05 

0.96766 

05 

0.96766 

05 

0.96766 05 

52 

L SC 1 1C 

0.5333E 05 

0.53336 05 

0.53336 05 

0.53336 05 

0.53,336 05 

0.53336 

05 

0.53336 

05 

0.53336 05 

53 

L SCI ID 

0.6578E 05 

0.65786 05 

0.6578 E 05 

0.6578E 05 

0.6578E 05 

0.65786 

05 

0. 65 786 

05 

0.S578E 05 

54 

LSC22T 

0.3714E 04 

0.37186 04 

0.3544E 04 

0.35646 04 

0.35666 04 

0.3572E 

04 

0. 35686 

04 

0.35766 04 

55 

ISC22C 

0.17066 05 

0.17886 05 

0.1704E 05 

0.17136 05 

0.1714E 05 

0.17176 

05 

0. 1 71 56 

05 

0.17196 05 

56 

LSC12 

0.2547E 04 

0,25476 04 

0.2547E 04 

0.2547E 04 

0.2547E 04 

0.25476 

04 

0.25476 

04 

0.25476 04 

57 

ISC 2 3 

0.18666 04 

0^ 1 866 6 04 

0.18656 04 

0. 1866E 04 

0.1 8666 04 

0.18666 

04 

0. 18666 

04 

0.18666 04 

58 

L SCC 23 

- 0. OOOOfc- 1 9 

-0.00006-19 

-0.3000E-19 

-0.00006-19 

-O.OCOOE- 19 

-0.00006 

-19 

-0.00006 

-19 

-0.00036-19 

59 

L SCC 13 

'0. 0003E-19 

-0.00006-19 

-0.03036-19 

-0 .00306- 19 

-0 .00006- 19 

-0.00006 

-19 

-0.00006 

-19 

-0.30006-19 

60 

LSCDF 

-O.OOOOE-19 

0.10526-01 

0.10526-01 

0.1052E-01 

0.10526-01 

0.10526 

-01 

0.10526 

-01 

0.1052E-01 

61 

Kl 12A6 

0.1371E 01 

0.13716 01 

0.1371E 01 

0 • 137 IE 01 

0.1371E 01 

0.13716 

01 

0.13716 

01 

0.13716 01 

62 

HUE i E 

-0. 23876 01 

-0.1400E 02 

-0.2529E 02 

-O.4504E 02 

-0 .68686 02 

-0. 10036 

03 

-0.12826 

03 

-0.17746 03 

63 

KfcLROT 

-O.OOOOE-19 

0.96056 00 

0.2757E 00 

-0.213 IE 00 

-0.95 18E 00 

-0.3178E 

01 

-0.94466 

01 

-0. 9204 E 01 

64 

EP S 11 

0.U32E-02 

0.29746-03 

3.2429E-02 

0 .278 IE-02 

0.4023E-02 

0.4375E 

-02 

0.65 06E 

-0? 

0.5572E-02 

65 

EPS22 

-0. 5B60E-02 

-0.67 086-02 

-0.1052E-01 

-0. 1256E-01 

-0.154BE-01 

-0.17526 
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(a) Fiber. Properties needed: Em 22, 33u (b) Matrix. Properties needed: E m u 22,33"' 

m2, 23, 13; G fl2, 23, 13: a fll, 22, 33i vml2, 23, 13; G fl2 23, 13; <>mll, 22, 33; 

Kfll, 22, 33; Hof, pf, N fi d fl - S f( . K m ll, 2233; H C m. Pm. Smc. e mpt e mpc . 

e mps' e mpT0R- 



tc> Ply. Input.- Fiber and matrix properties: Pg, p^, ps, 
AT. Properties computed: Em_ 22, 33; VU2, 23 13; 
g U 223,13; am, 22 33; Kai.22,33; H c l H <1 
s aiT, 11C, 22T, 22C, 12S, 23S; *112. Stress analysis: 
e Ul,2212 °U1, 22 12> 1.0 - F(o,S,K n2 ). 


Ply 

orientation 

angle, 

9(K deg z 



Figure 1. - Typical multilayered fiber composite and some basic definitions. 
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Figure 2 . - Code physical arrangement. 




Figure 3. - Physical arrangement of input data cards, 


















E-glass /epoxy 





































Figu re 6. - Ply orientation geometry. Composite structural axes, x,y,Z; composite 
material axes, 1,2,3; ply material axes (coincides with fiber direction, lj, 2j, 3j). 
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